Radiation of Terahertz Electromagnetic Waves from Build-in Nano Josephson Junctions of Cuprate 

High-r c Superconductors 
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The nano-scale intrinsic Josephson junctions in highly anisotropic cuprate superconductors have potential for 
generation of terahertz electromagnetic waves. When the thickness of a superconductor sample is much smaller 
than the wavelength of electromagnetic waves in vacuum, the superconductor renders itself as a cavity. Unlike 
conventional lasers, the presence of the cavity does not guarantee a coherent emission because of the internal 
degree of freedom of the superconductivity phase in long junctions. We study the excitation of terahertz wave by 
solitons in a stack of intrinsic Josephson junctions, especially for relatively short junctions. Coherent emission 
requires a rectangular configuration of solitons. However such a configuration is unstable against weak fluctua- 
tions, contrarily solitons favor a triangular lattice corresponding to an out-phase oscillation of electromagnetic 
waves. To utilize the cavity, we propose to use an array of stacks of short intrinsic Josephson junctions to gener- 
ate powerful terahertz electromagnetic waves. The cavity synchronizes the plasma oscillation in different stacks 
and the emission intensity is predicted to be proportional to the number of stacks squared. 
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INTRODUCTION 

Generations of electromagnetic (EM) waves are mainly 
based on two principles; the low frequency waves are gener- 
ated by electronic oscillations while the high frequency waves 
are generated by photonic methods. Because of the fundamen- 
tal electron-velocity limits, the performance of microwave 
electronic circuits degrades very rapidly above 100GHz. The 
photonic generators are also limited to the frequency higher 
than several terahertz (THz). There exists a frequency gap 
ranging roughly from O.lTHz to lOTHz, which is usually 
dubbed as the THz gap. The THz waves, on the other hand, 
have potential for wide applications, such as drug detection, 
security check and so on|l|,|2|]. Thus it is extremely demand- 
ing to develop solid-state, compact generators of THz EM 
waves. The Josephson effect in cuprate high-r e superconduc- 
tors is promising for this purposeyj. 

When two superconductors get sufficiently close to each 
other to form a tunnel junction or Josephson junction, macro- 
scopic quantum phase coherence will be established between 
them. According to the uncertainty relation between the num- 
ber of Cooper pairs and the superconductivity phase, a super- 
current appears in the tunnel junction, which is given by the 
dc Josephson relation^ 



/ = LsinP, 



(1) 



with P being the phase difference between these two super- 
conductors and I c the critical current. Biased by a dc voltage 
V, the Josephson junction is essentially a two-level system 
with energy difference 2eV. As Cooper pairs are tunneling, 
EM wave is emitted from the junction, whose frequency is 
governed by the ac Josephson relation 



hco = hd,P = 2eV. 



(2) 



For instance, lmV corresponds to 0.483THz. In principle, the 
operating frequency can be tuned continuously by the biased 



voltage. In contrast to conventional lasers, the population in- 
version can be achieved in a meta-stable superconducting state 
with finite normal resistance realized in a junction with small 
conductance. The radiation from a single junction, however, 
is only about lpW and the frequency is below THz limited by 
the small superconducting energy gapQQ]. 

One may fabricate an array of Josephson junctions to en- 
hance the radiation power||3 Actually, it was found 
in 1992 that layered high-T c superconductors, such as 
Bi2Sr2CaCu208 + a(BSCCO), intrinsically form a stack of 
Josephson junctions (IJJs) of nano-scale thickness^, as 
drawn schematically in Fig. \V[a). The advantages of IJJs 
are: 1) large superconductivity energy gap, which allows 
for operations up to 15THz; 2) almost homogeneous in the 
atomic scale, which makes superradiation possible; 3) these 
nano-scale build-injunctions render a huge inductive coupling 
between neighboring junctions. The proposed device using 
BSCCO single crystal to generate THz wave is depicted in 
Kg. Hlb). 

The following two conditions must be met in order to 
achieve strong THz emission. First, coherent emission of EM 
waves involves synchronized oscillations of individual ele- 
ments. Synchronization requires mutual couplings between 
each oscillator and usually involves long-range interactions. 
In the context of the Kuramoto model lllOll . it was found that 
the lower critical dimensions for long-range synchronous os- 
cillations is four ll ill . Practically, one introduces a common 
resonator that couples to all oscillators. This resonator pro- 
vides a feedback mechanism and establishes a long-range in- 
teraction between oscillators, which forces all oscillators to 
run at its resonant frequency. For one dimensional array of 
Josephson junctions connected to a common RLC resonant 
circuit, it was demonstrated rigorously that the system can 
be mapped to the infinite dimensional Kuramoto model 11211 . 
where each junction is coupled to the rest junctions. For a 
stack of IJJs, it was formulated theoretically OGIl] that there 
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exists a significant impedance mismatch for a stack of thick- 
ness of several micrometers as in the recent experiments lfl5l 
[lcjll . similar to that in a thin capacitor. Thus a thin stack of IJJs 
intrinsically forms a cavity, which is helpful for synchronizing 
all junctions. 

Secondly, for a stack of IJJs to work as an efficient gain 
medium, the supercurrent should be modulated in the lat- 
eral directions such that considerable dc input power can be 
pumped into the Josephson plasma oscillation, a collective 
composite wave of electromagnetic oscillations and back-and- 
forth tunneling of Cooper pairs. One well-known example of 
the modulation in a single Josephson junction is 2n change 
in the superconductivity phase, corresponding to a shuttling 
soliton, a nonlinear quantized topological excitation. An- 
other example is the recently discovered new dynamic state 
in a stack of strong inductively coupled IJJs ll 171 - 12 ill , where n 
phase changes are localized at the nodes of electric fields and 
stacked periodically along the c axis. By exciting the in-phase 
plasma mode, this n kink state supports strong THz emissions. 

Yet a cavity becomes a standard ingredient in most lasers 
and is almost indispensable to coherent radiation. However, 
for a stack of long IJJs embedded in the cavity formed by the 
stack itself, because of the existence of internal degree of free- 
dom, the superconductivity phases of different junctions may 
organize themselves in a way that the oscillations in neighbor- 
ing junctions are out-phase. To be concrete, let us consider the 
stack shown in Fig. [TJ a )- The stack forms a two dimensional 
cavity and the plasma mode is cosimjix/ L) smiqnzl L z ) or is 
simply denoted as (m, q) mode. Here L is the length of the 
stack and L z is the thickness. Coherent emission from edges 
at x = or x — L requires the excitation of q — 1 mode which 
is not realized in some specific configurations of the super- 
conductivity phase such as solitons as will be discussed later. 
Therefore it is of both theoretical and practical importance to 
understand the behavior of superconductivity phase in a reso- 
nant cavity. 

In the present paper, we investigate the behaviors of solitons 
in a stack of IJJs, especially in a relatively short one. We find 
that in a short stack with radiating boundary condition, the 
soliton becomes incomplete and resembles oscillating stand- 
ing waves. We also find that the rectangular lattice of solitons 
associated with in-phase plasma oscillation is unstable. The 
out-phase oscillation of soliton in a cavity leads us to consider 
an array of short junctions enclosed in a cavity to avoid the 
internal structure of phases in individual stacks. This config- 
uration provides an alternative scheme to achieve strong THz 
radiation from Josephson junctions. 



MODEL 

We assume that the stack is infinitely long in the y direction 
and the system is reduced into two dimensions as shown in 
Fig. QI a )- The dynamics in a stack of IJJs is governed by the 
well-known inductively coupled sine-Gordon equationsj22j- 
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FIG. 1: (a) Highly anisotropic cuprate superconductors like BSCCO 
can be modeled as a stack of IJJs. The blue (gray) region denotes 
superconducting (insulating) layers, (b) Proposed device to generate 
THz EM waves. The BSCCO single crystal is sandwiched by two 
gold electrodes through which a dc current is fed in. 
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d 2 x P, = (1 - <rA d )[sinP, + B3.P, + d)P x - /ex,], (3) 

where P\ is the gauge invariant phase difference at the Z-th 
junction, £ the inductive coupling of order of 10 5 , B the nor- 
malized conductivity and J ext the external current i20ll . Length 
is in unit of the London penetration depth A c ss 200/im and 
time in unit of the Josephson plasma frequency u>j » 0.5THz. 
Ad is the finite difference operator defined as Ad/i = fi+i + 
/ M - 2/,. We use B = 0.02 and f = 7.1 x 10 4 , which 
are typical for BSCCO. I20I1 The normalized magnetic field 
Bi and electric field E\ in Z-th junction are related to Pi by 
(1 - £Ad)£) = d x Pi and E\ = d t Pi. Because of the signif- 
icant impedance mismatch, the oscillating magnetic field is 
extremely small at the edges, thus as a good approximation, 
the boundary condition can be written as 

d,P, = Zd x P h (4) 

with the impedance of the stack Z ~ Aem/L z » 1. Here 
^em ~ 300yum is the wavelength of THz EM wave in vacuum 
and L : ~ l/vm is the thickness of the stack. We consider a 
situation where no external magnetic field is applied to the 
stack. Equation (f3]i together with the boundary condition Eq. 
(0 is solved numerically and the electromagnetic fields are 
then derived from P/. 
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FIG. 2: (a) Time evolution of P; in a long junction L = 4A C and with 
the closed boundary condition Z = oo. (b) Time evolution of Pi in 
a relatively short junction L = 2A C and with radiating boundary con- 
dition Z = 10. The soliton under this situation becomes incomplete, 
(c) Frequency spectrum of the electric field at a voltage of V = 2.79 
associated with the first zero-field step and 7 cxt = 0.1 for L = 2A C . 
The fundamental peak does not come to the frequency predicted by 
the ac Josephson relation from the voltage. 

SOLITONS SOLUTION AND ITS STABILITY 

Equation ([3]) exhibits rich dynamical behaviors^ and one 
possible dynamic state is the solitons state. A soliton car- 
ries quantum flux <I>o = hc/2e as 2.0 x 10~ 15 Wb, and it ex- 
periences repulsive force from other solitons with the same 
polarity. The interplay between the interaction between soli- 



tons and the cavity gives rise to rich configurations of solitons, 
from the rectangular lattice to triangular one. For a rectangular 
configuration of solitons along the z axis, Eq. (0 is decoupled 
and the problem essentially becomes single-junction physics. 
We first discuss the properties of the uniform soliton solution 
and then consider its stability. 

A soliton has characteristic length of A c and undergoes 
Lorentz contraction when its velocity approaches the plasma 
velocity, i.e. the light velocity in junctions. Thus, to support 
a pronounced soliton structure, the length of junction should 
be much larger than A c . As shown in Fig. |2ja) for L = 4A C , 
a well-defined soliton propagates leftwards and get reflected 
at the left edge and becomes an anti-soliton. This anti-soliton 
then travels rightwards and repeats the above process. In one 
period, the phase increases by 2n and the period is given by 
the voltage according to Eq. (|2]). The periodic motion and re- 
flection of solitons inside the cavity excite plasma oscillation 
at a frequency depending on the pattern of solitons in the x di- 
rection, and induce current steps at voltages V = 2nn/L with 
n being the number of solitons in the lateral direction [20]. 

It would be interesting to study the structure of soliton in a 
junction with length comparable to A c and with slightly open 
boundary condition. In this case, it is difficult to see a com- 
plete soliton traveling across the junction and the dynamics 
of Pi resemble oscillations of standing waves in addition to a 
linear term in time, as shown in Fig. |2{b). However, because 
of the commonly underlying parametric instability for soli- 
tons, the frequency spectrum is the same as that of long junc- 
tion, namely many frequency harmonics are visible, and the 
radiation frequency and voltage does not obey the ac Joseph- 
son relation as depicted in Fig. [2jc), opposed to the recent 
experiments iflU . 

In a recent work, a new state exhibiting standing wave os- 
cillations but with the same IV characteristics as those of soli- 
tons solution was found in a relatively short junction J2g]. Our 
simulations of solitons in a junction with length comparable 
to A c reproduce all the results in Ref. If26ll . which lead one to 
assign the dynamic state in Ref. If26ll as a soliton state. 

We then investigate the stability of the uniform soliton so- 
lution by including small fluctuations into the system, e.g. 
adding white-noise currents at the right-hand side of Eq. ([3]). 
The uniform solution then becomes unstable and the system 
evolves into the n kink state. The uniform solution cannot 
be stabilized by introducing dissipation, radiation or a small 
magnetic field either. Indeed, due to the strong mutual re- 
pulsive interaction, the solitons favor triangular lattice corre- 
sponding to the q — N cavity mode with N the number of 
junctions in the stack. Thus in the soliton state, due to the 
internal degree of freedom, the cavity does not guarantee a 
coherent oscillations in all junctions. 

AN ARRAY OF STACKS OF IJJS 

To utilize the cavity, one way is to adopt an array of stacks 
of short IJJs to avoid the complicated behaviors of Pi in indi- 
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vidual stacks. A schematic view of such a setup is depicted in 
Fig. Ha). 

We consider short IJJs with L <K A c , thus the soliton state 
cannot be stabilized. The short IJJs also allow for efficient 
cooling and less Joule heating [27]. For simplicity, the height 
of the cavity is assumed to the same as of IJJs and L v is as- 
sumed to be infinite as before. For L z <s; A^m, only the mode 
with the wave vector k z — can propagate in vacuum. This 
uniform mode tends to excite in-phase mode in IJJs. The elec- 
tromagnetic waves in the cavity are therefore governed by 

d 2 x P = J' c {x) sin P + p'(x)d,P + e'(x)d 1 P - J' en (x). (5) 

In IJJs, Eq. <(5j is the uniform-mode correspondence of Eq. 
©, while in vacuum P should be interpreted as an auxiliary 
variable determined by d x P - B y and d,P - E z . Injunctions 
p'(x) = p, J' c (x) = J c , J' e Jx) = 7 ext and e'{x) = 1, while 
in vacuum yS'(x), J' c (x), J' ext (x) are zero and e'(x) = ej with e,/ 
being the dielectric constant of vacuum. We assume a metallic 
cavity, and at the edges of the cavity x = or x — L c , we 
have E z = for an ideal cavity. The radiation effect is taken 
into account by introducing an effective surface impedance Z c 
for the cavity. Then the radiation power is given by S r = 
\E 2 Z \/(2Z C ). 

The solution to Eq. (0 in the linear region is 

P - cot + Re[-iAg(iL> c , x) exp(/w?)], (6) 

where the first term is the rotating phase responsible for the 
dc voltage and the second term is the uniform plasma oscil- 
lations along the z axis. g(u> c , x) is the cavity mode with the 
resonant frequency 0) c , which can be determined by neglect- 
ing the Josephson current, external current and dissipations, 
since the Josephson current and external current are responsi- 
ble for the excitation of EM waves while dissipation has the 
effect of broadening the resonance line width. The profile of 
g(x) for the first mode is shown in Fig. Ob), where there is a 
half wavelength residing in the cavity. The amplitude of the 
plasma oscillation A then is given by 

A = Iq C dxJ' c (x)g(x) 

jf c dx [(ofi - a%) e'(x) - i/3'(x)<o] g(x) 2 ' 

The radiation slightly modulates the supercurrent in each 
stack and pumps energy into the plasma oscillation. The 
amplification of the plasma oscillation is most efficient for 
the first cavity mode because the radiation in the region with 
g(x) < is anti-phase for higher modes. For the first mode, 
we have A ~ Nj with Nj being the number of IJJs stacks. Thus 
the total radiation power increases as Nj squared as shown in 
Fig. [3tc), which presumes a superradiation. It is remarked 
that IJJs are not necessarily identical since the cavity helps to 
synchronize them. 10,0] 
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FIG. 3: (a) Schematic view of an array of IJJs embedded in an ex- 
ternal cavity, (b) Profile of g(x) for the first cavity mode, (c) De- 
pendence of the radiation power S r ~ \cl>A\ 2 /(2Z C ) on the number of 
IJJs stacks Nj. In obtaining the results in (b) and (c), we assume that 
IJJs stacks with length 10//m are evenly distributed in the cavity with 
dielectric constant €d =0.1. Here L c = 400/rni, and Nj = 8 in (b). 



conductors and the associated dynamics of superconductivity 
phase. We studied the behaviors of a soliton in a relatively 
short junction, and found that the soliton becomes incomplete 
and resembles the oscillating standing waves. We also found 
the in-phase configuration of solitons is unstable even though 
they are embedded in a cavity. The instability is caused by 
the mutual interaction between solitons. Based on these ob- 
servations, we proposed to integrate an array of stacks of short 
intrinsic Josephson junctions in an external cavity to attain the 
powerful terahertz emission. The short junctions avoid com- 
plicated behaviors of the superconductivity phase and allow 
for efficient synchronization by the external cavity. Once they 
are synchronized, the radiation power increases as the number 
of stacks squared according to our calculations. 
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CONCLUSIONS 



To summarize, we have discussed the possible radiation of 
terahertz electromagnetic waves from cuprate high-T c super- 
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